Improvement via hypercubic smearing in triplet and sextet QCD 



Yigal Shamir, Benjamin Svetitsky, and Evgeny Yurkovsky 
Raymond and Beverly Sackler School of Physics and Astronomy, Tel Aviv University, 69978 Tel Aviv, Israel 

We study non-perturbative improvement in SU(3) lattice gauge theory coupled to fermions in the 
fundamental and two-index symmetric representations. Our lattice action is defined with hypercubic 
smeared links incorporated into the Wilson-clover fermion kernel. Using standard Schrodinger- 
functional techniques we estimate the clover coefficient csw and find that discretization errors are 
much smaller than in thin-link theories. 
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I. INTRODUCTION 

The improvement of a lattice action is meant to reduce 
the effects of lattice artifacts and thus to bring calculated 
quantities closer to their continuum limits. In our work 
on the SU(3) gauge theory with sextet fermions [H-ll] we 
have adopted normalized hypercubic (nHYP) smearing 

^] in the expectation that it would yield reliable re- 
sults on fairly coarse lattices. We found, indeed Q, that 
the smeared-link action keeps the critical hopping pa- 
rameter Kc(/3) of the Wilson-clover action much closer to 
its continuum value of 1/8 even for strong bare couplings 

= 6//3. It also pushes a first-order phase transition to- 
wards stronger bare couplings, thus allowing calculation 
of the running coupling g^y in a regime where both 
and (7gp are quite strong. Moreover, the smeared plaque- 
tte averages are much closer to unity than the "thin-link" 
plaquettes, a clear sign that the gauge field is being ef- 
fectively smoothed. 

In this note we show that nHYP improvement extends 
to the axial Ward identity (AWI). In a calculation in- 
volving Wilson fermions, the AWI is frequently used to 
determine Kc by demanding that the quark mass m be 
zero, a consequence of conservation of the isovector ax- 
ial current. One can adjust the improvement coefficient 
Csw, which multiplies the clover term in the fermion ac- 
tion [7|, in order to minimize errors in the AWI itself. 
One sign of such errors is the sensitivity of m to the 
location where the AWI is measured. Liischer and col- 
laborators Q proposed a procedure of measuring m at 
points of the lattice that are inequivalent because of the 
boundary conditions used in the Schrodinger Functional 
(SF) method. Demanding that m = at two such points 
gives conditions for fixing Kc and csw- 

Liischer et al. applied the AWI criterion to the 
quenched SU(3) gauge theory, defining the axial current 
using thin links. While csw = 1 at tree level, they found 
numerically that cgw — 1-8 is required for = 6. Jansen 
and Sommcr [10] did a similar calculation for QCD with 
thin-link Wilson fermions {Nf — 2) and again found that 
a large value of csw — 1 is required. Defining the current 
with nHYP-smeared links, Hoffmann, Hasenfratz, and 
Schaefer [ll| reduced the required csw in the quenched 
theory from 1.55 to 1.05 at /? ~ 6.4. This is dramatic 



evidence for the claim that nHYP smearing brings the 
theory much closer to the continuum limit. 

In this paper we present a calculation in theories with 
Nf — 2 dynamical triplet and sextet quarks. For quarks 
in the triplet representation, we compare the thin-link 
and smeared-link theories and find results as dramatic 
as in the quenched theory. We find small violations of 
the AWI in the smeared-link sextet theory as well. We 
conclude that setting csw to 1 is adequate when nHYP 
smearing is used. 

The procedure of keeping csw — 1 while using smeared 
links was advocated for triplet QCD in Ref. [l3| and 
subsequently tested in large-scale calculations |13l . IT3 |. 
For a determination of csw with partial stout smearing, 
see [H. 



II. DETERMINATION OF csw 

We follow closely the jiethod of Refs. The 
0(a)-improved action is 



S ^ S\N[U,^,-ip]+ a^csw^i^^u^^F^^i), (1) 



where S'w is the conventional Wilson action comprised of 
the plaquette gauge action and the fermion hopping term. 
The second term in Eq. ([T]) is the clover term, wherein 
F^^ is the lattice-discretized field strength. The gauge 
links in the hopping and clover terms are "fat links," 
defined via nHYP smearing as described in Ref. @ and 
with the same smearing parameters. 

On a lattice with x T sites, we impose Dirichlet 
boundary conditions on the gauge fields on the temporal 
boundaries of the lattice. 



2 



with the asymmetric choice"'^ 

Ck = — diag(-7r, 0, tt), = — diag(-57r, 27r, Stt). 

(3) 

The fermion boundary conditions are homogeneous, 





P_V'(i = 0) 


= p+V(t = r) =0, 






^{t = 0)P+ 


= Tjj{t = T)P^ =0, 


(4) 


with Pi - 


- 2 V ' • 


The gauge-invariant boundary 


fields that 


can be used 


as wah sources are 




c - 


X 


0)P+ip{x,t^ 1) 


(5) 


c - 


X 


= l)P_[/t(x,t = 0) 


(6) 


c - 


X 


= T-l)P_?/;(x,i = T-1) 


(7) 




X 


= r-l)P+C/o(x,t = T-l). 


(8) 



(These are the same as those used by 8], but written 
in exphcit form.) The gauge fields are periodic in the 
spatial directions, while the fermion fields satisfy ip{L) ~ 
e"/5?A(0) [81 . 

The axial Ward identity 



d,AZ,^2mP'^ 



(9) 



gives a definition of the quark mass m. Equation Q 
contains the 0(a)-improved axial current. 



(10) 



The pseudoscalar density and the unimproved axial cur- 
rent are defined by the local products 

p- = yV, a; = ^n^ij-^i^. (11) 

In practice, one evaluates correlation functions of Eq. ^ , 
d,{^:,^{x)d) = 2m{P%x)d), (12) 

where O is any operator located at non-zero distance 
from X. It is convenient [1] to use the pseudoscalar field 
made of the boundary operators to define wall sources O 
and O' at i = and t — T, respectively, viz., 



0" = C75yC, O'-'-C'TsyC'. (13) 



Note the appearance of T in Eq. The authors of Refs. 

use L instead, which means that for the lattice used here our 

background field is weaker than theirs. 



Thus we define the correlation functions 

1 

3^ 
1 

3^ 



fpixo) 

fA{xo) 



{P''ixo)0-) 



(14) 
(15) 



which depend only on xq by translation invariance of the 
wall sources. The spatial derivatives in Eq. (|12p vanish 
similarly, whence one obtains an estimate for m. 



m{xo) 



dpfAjxp) + CAadofp{xo) 
2fp{xo) 



(16) 



where for do we use a symmetric derivative and dg is the 
nearest-neighbor second derivative. Using O' we define 
analogously 



fp{T-xo) = -l{P''{xo)d'^) 



}'a{T - xo) 



--(Ag(xo)O"^), 



(17) 
(18) 



which leads to an alternative estimate m'(a;o), defined by 
the parallel of Eq. (HH) in terms of f'p, /'a- If the AWI is 
respected by the improved action then one would expect 
that m and m' are independent of xq and equal to each 
other. 

Equation (|16p and its primed counterpart still contain 
the unknown coefficient ca- An alternative definition of 
the mass eliminates this dependence. If we write Eq. 
as 

m{xo) = r{xo) + CAs{xo), (19) 

and similarly for m' ,r' , s' , the alternative is 

,.r/ X / N / -.riyo) -r'{yo) , . 

M{xo,yo) = r{xo) - s{xo)-— — — , (20) 

s(yo) - s'(yo) 

which differs from m,m' in O(a^). We also define the 
quantity M' by exchanging primed and unprimed vari- 
ables in Eq. Then, still following Refs. [i-[ll|, we 
define the measure of residual violation of the AWI to be 



AM = M ( ^T, -^t\ - M' (-^T, 



III. NUMERICAL RESULTS 



(21) 



All our calculations were performed on lattices with 
i = 8, T = 16. For triplet as well as for sextet quarks 
we chose two values of /3, in the weak and intermediate 
coupling regions. Varying cswi we calculated AM as 
defined in Eq. ((2T|) . Our results are displayed in Tables H]- 
HVland in Figs. [Hand H 

For the fat-link theories, we fixed k for each f3 by de- 
manding that r{T/2) = at cgw — 1, that is, by setting 
to zero the unimproved quark mass. This is an alterna- 
tive to requiring, say, M{xo, j/o) = for some Xq, yo; it fol- 
lows on the observation that the second term in Eq. ([20)1 
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/3 


K 


csw 


aAM X lO" 


7.4 


0.1255 


1 


1.3(11) 






1.1 


-3.2(7) 






1.2 


—9.3(14) 






1 Q 
i.O 


— io.o^ioj 


5.8 


0.1267 


1 


7.2(73) 






1.1 


-1.1(23) 






1.2 


-6.6(19) 






1.3 


-16.4(31) 



TABLE I: Results of AM calculations on 8 x 16 lattice. 
SU(3) gauge theory with Nf = 2 fermions in the fundamental 
representation, fat links. 2000 trajectories were run at /? = 7.4 
for each value of csw, and 3000 at /3 = 5.8. 
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-5.4(22) 

39.(15) 

8.5(30) 
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TABLE II: AM for thin-link fermions in the fundamental 
representation, for comparison with Tabled Values of /3, k, 
and Csw were chosen as in Ref. [lo| . 4000 trajectories were 
run at P — 7.4 for each value of csw, and 18000 at /3 = 5.7. 



is generally small and recognizes the fact that the AWI 
([9|) should hold for nonzero mass as well. In comparing 
to other work we note that the K-dependence of AM has 
been seen to be very weak 0, . For the same reason, 
we did not vary k for the fat-link theories as we varied 
Csw at given /3. 

The optimal value Cgw(/3) is determined by demanding 
AM = 0. We do this via linear fits to AA/(csw), with 
the results shown in Table ITVl and plotted in the figures. 

In the theory with triplet fermions we compare the 
fat-link results to AM in the thin-link theory, calculated 
at values of /3, k, and csw used by Jansen and Som- 
mer [lo| . [In this case, k was shifted with csw to keep 
M(iT, \T) = 0.] It is clear that for csw - 1 the fat-link 
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TABLE III: As in TableU but for fat-link fermions in the sex- 
tet representation. 4000 trajectories were run for each /3, csw- 



FIG. 1: Determination of csw for the fat- and thin-link theo- 
ries with triplet quarks, using the AM data in Tables HI andUll 
The curves are linear fits whose slopes and intercepts are given 
in Table HVl 




FIG. 2: Data for the fat-link theory with sextet quarks (Ta- 
ble |III|), plotted together with the fat-link triplet data from 
Fig. [1] The linear fits shown have the parameters listed in 
Table IIVI Data points for [3 = 5.8 have been given slight 
horizontal shifts for clarity. 



action gives much smaller values of AM for both values 
of /?. The fat-link action also gives Cgw(/3) very close to 
unity, even at the stronger coupling (3 = 5.8, for both the 
triplet and the sextet theory. 



IV. DISCUSSION 

Figs. [T] and O and Table IIVI show that discretization 
errors in both fat-link theories are generally small, as re- 
flected in the value of AM at csw — 0. In comparing 
thin- and fat-link fermions for the triplet theory, we find 
that the slope of AM vs. csw is similar but that the non- 
perturbatively determined coefficient is much closer 
to one for the smeared links. 
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-0.005(1) 
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-0.0057(5) 


/3 = 


5.7, triplet 
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5.8, triplet 
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-0.007(2) 
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1.07(2) 


0.0033(5) 


/3 = 


5.8, sextet 


fat links 


1.17(6) 


0.002(1) 



TABLE IV: Results of fitting AM(csw) to a straight line for 
each theory: the value of csw where AAf = 0, and the fitted 
slope. 

Our comparison between thin- and fat-link theories 
uses results obtained at (roughly) the same values of the 
bare coupling. One could compare instead at points of 
similar physics by demanding that the values of /3 give, 
for example, the same string tension. This would consti- 
tute a comparison of the different discretizations at equal 
lattice spacing. Our results, however, make this elabo- 
rate exercise unnecessary. In contrast with the triplet 
thin-link results, the fat-link results change very little in 
going from /? = 7.4 to 5.8. Thus, had we tuned our bare 
couplings to the slightly different values needed to repro- 
duce the lattice spacings of Ref. [10|] , there would be no 
qualitative change in the conclusions. 

Turning to the sextet theory, one might be concerned 
by the small slope of AAf(csw), which acts to enlarge 
the error bar on Cg^. Is this due to having chosen a 
quantity with low sensitivity to csw? Let us consider 
the tree-level value of AAf, denoted AM^^^. This quan- 
tity sets a scale for AWI violation by the background field 
on the finite lattice. It is easily determined numerically, 
by calculating the Green functions in Sec. In] in the single 
gauge configuration that is the classical solution to the 
boundary conditions, namely, a constant electric field. 
(We set Csw = 1 and k = 1/8.) In the triplet theory we 
obtain aAAf(°' = 0.00025, whereas for the sextet theory, 

^ For a similar determination of cgw in the SU(2) gauge theory 
with thin-link fermions in the fundamental and adjoint represen- 



aAAf = -0.00065. In the sextet theory the fluctu- 
ations of the dynamical fields weaken AAf at csw = 1 
considerably from its tree-level value. The observation 
that AAf (csw = 1) is small on the scale set by AAf (''^ 
means that there is no reason to vary csw away from 1. 

The authors of Refs. 0-[IH fixed by demanding 
AAf = AAf(°^. This is an attempt to determine Cs-^ 
in infinite volume, by supposing that the finite-volume 
corrections to the dynamical AAf and to the tree-level 
AAf are the same. Since this assumption is not yet 
supported by any study of volume dependence, we eschew 
this procedure in favor of demanding AAf = 0. In other 
words, we determine Cg^ for our volume without making 
any statement about the L — )• oo limit. 

The smallness of discretization effects we report for the 
fat-link theories is consistent with the reduced discretiza- 
tion errors found in other observables in the SU(3) /sextet 
theory P as well as in the SU(2)/adjoint theory [l6|.^ 
This lends support to our decision to stick with the tree- 
level value, Csw = 1- Indeed, thanks to asymptotic free- 
dom, the continuum limit of a lattice gauge theory with a 
given (Dirac-) fermion content is completely determined 
once the fermion masses are fixed. Improvement is the 
art of reducing discretization errors when the lattice spac- 
ing is nonzero. Since it does not change the continuum 
limit, improvement is always at one's discretion, never 
mandatory. Our experience is that, thanks to the fat-link 
action, using csw = 1 provides as much improvement as 
we need. 
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